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Abstract 
o 

c/3 ' We review recent ah initio molecular dynamics studies of electrode/electrolyte interfaces in 

lithium ion batteries. Our goals are to introduce experimentalists to simulation techniques ap- 
plicable to models which are arguably most faithful to experimental conditions so far, and to 
rH \ emphasize to theorists that the inherently interdisciplinary nature of this subject requires bridging 

F^ ' the gap between solid and liquid state perspectives. We consider liquid ethylene carbonate (EC) 

Q . decomposition on lithium intercalated graphite, lithium metal, oxide-coated graphite, and spinel 

manganese oxide surfaces. These calculations are put in the context of more widely studied water- 
t^ solid interfaces. Our main themes include kinetically controlled two-electron-induced reactions, the 

r~^ ■ breaking of a previously much neglected chemical bond in EC, and electron tunneling. Future work 

o\'. 

if^ , on modeling batteries at atomic lengthscales requires capabilities beyond state-of-the-art, which 

^s I emphasizes that applied battery research can and should drive fundamental science development. 

m 
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I. INTRODUCTION: BACKGROUND AND CONTEXT 

Lithium ion batteries (LIB) are currently the devices being implemented or considered 
for large scale static and transportation energy storage. They carry high energy density and 
have the potential of dramatically reducing green house gas emission because they operate 
within a high voltage window.^ Today's commercial LIBs (Fig. [1^) consist of graphitic car- 
bon anodes, transition metal oxide cathodes, and organic solvent-based electrolyte. Other 
crucial LIB components include passivating "solid electrolyte interphase" (SEI) films formed 
from excess electron-induced electrolyte decomposition products on anode surfaces.-"- SEI 
films are heterogeneous in structure and consist of Li2C03, ethylene dicarbonate (EDC), 
oligomeric/polymeric compounds, salt decomposition fragments, and other products.-"- 
They prevent continuous electron injection into the electrolyte, averting further loss of Li"*" 
and electrolyte molecules. Li"*" transport through SEI films remains adequately fast. Oxi- 
dation products are also often found on cathode surfaces. New concepts of electrodes being 
pursued, such as Si-based anodes and "air" cathodes in metal-air batteries, share many 
solid-liquid interface features shown in Fig. [1^. LIBs are pragmatic devices. They combine 
our best electrochemical, solid, and liquid state expertise to deliver high volumetric and 
gravimetric energy/power densities. While all-solid batteries have received much attention 
for niche applications and all-liquid flow-cell batteries have significant potential for static 
storage, batteries featuring both liquid and solid components will undoubtedly dominate for 
the foreseeable future. 

It has been widely acknowledged that interfaces are critical for good performance and 
long lifetime in batteries.- To some extent, interfaces dictate the choice of electrode ma- 
terials and electrolytes. Graphitic anode forms stable SEI with ethylene carbonate (EC), 
not propylene carbonate.'^ Another celebrated example of the interconnected nature of LIB 
degradation concerns spinel lithium manganese oxide cathodes. Mn(II) ions dissolve from 
the spinel^ and diffuse to the anode. They are incorporated into the SEI there and degrade 
its passivation properties via mechanisms not yet fully understood, leading to enhanced re- 
ductive decomposition of the electrolyte.- Thus solid and liquid degradation modes can be 
strongly coupled. The choice of carbon materials also strongly influences their viability as 
anodes.'* This is possibly related to interfacial effects such as instability towards exfoliation 
of graphite induced by solvent intercalation^ and carbon-edge functional groups. Molecular 



additives like vinylene carbonate (VC) have been added to improve electrode passivation,- 
and these can react on both cathode and anode surfaces.- Interfaces are particularly per- 
tinent to nanostructured electrodes for energy storage applications where the large surface 
areas call for enhanced stabilization.— i^i If adequate control of interfaces can be attained, 
nanostructures have the potential of achieving much higher electron and lithium ion trans- 
port rates, and alleviating strain-induced electrode cracking and degradation. 

To some extent, theory and experiments used to study LIB liquid-solid interfaces are 
complementary. Experimental techniques that have been applied to SEI studies include 
Fourier Transform infrared spectroscopy,— electrochemical impedance spectroscopy,— X-ray 
photoelectron spectroscopy (XPS),^^"^^ atomic force microscopy and Raman spectroscopy,— 
nuclear magnetic resonance,— and transmission electron microscopy.— This article only sam- 
ples a few experimental papers and focuses instead on recent atomic-lengthscale modeling 
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of electrochemical reactions at LIB electrode/electrolyte interfaces.—"— (See Ref. 
comprehensive review of experimental methods and results up to 2004.) Modeling reactions 
involve predicting electron transfer and chemical changes and typically relies on Density 
Functional Theory (DFT) based techniques. Such calculations probe smaller length and 
time scales than experiments, and can reveal insights about thermodynamics and kinetics 
which are difficult to measure under the non-equilibrium, kinetically-driven SEI formation 
conditions. Theoretical SEI growth mechanistic studies were pioneered by Balbuena et al.^s 



cluster-based work^^ (see Ref. 
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for a theoretical overview). With the continual growth of 
computational power, the admittedly costly DFT-based ab initio molecular dynamic method 
(AIMD, discussed below) will be a mainstay of future liquid-solid interface studies. AIMD 
permits inclusion of electrodes and liquid electrolytes in model systems and adds multi- 
electron reactions, reactive surface sites, and new perspectives to Balbuena et a/.'s work. 
Other interfacial phenomena, like the dynamics of Li"'' intercalation between electrolyte and 
electrodes, are important but arguably must await elucidation of how electrode surfaces are 
modified by electrolyte decomposition products, unless the electrode operates at sufficiently 
modest voltages to keep its surfaces pristine.— i^^ 

Batteries are complex systems. The applied voltage, state-of-charge (lithium content), 
temperature, electrode cracking, sweep rate, presence of conductive carbon materials, native 
surface films like lithium carbonates on cathodes, and interference from salts/contaminants 
all affect interfacial behavior. They are challenging to model even on bare oxide surfaces 




FIG. 1: (a) A simple schematic (not an actual calculation) of a lithium ion battery depicting 
a graphitic carbon anode, a LiMn204 cathode, and ethylene carbonate (EC)-based electrolyte. 
During charging, Li^ moves from the cathode to the anode, accompanied by electron flow. A 
passivating SEI film, as labeled, starts to grow upon charging at the anode-electrolyte interface. 
Electrolyte oxidation products also emerge on the cathode surface, (b) EC with atomic labels, (c) 
Dimethyl carbonate (DMC). Grey, red, white, blue, and purple spheres represent C, O, H, Li, and 
Mn atoms respectively. Liquid state EC are depicted as stick figures. 

(i.e., in vacuum).—"— We have adopted a basic science approach. Our starting points are 
defect-free surfaces that can be probed using AIMD methods at finite temperature. Far from 
being routine DFT/AIMD apphcations, we will show that modeling clean battery interfaces 
already requires going beyond state-of-the-art theoretical capability and is a good example 
of application-driven fundamental research. 

Modeling liquid/solid interfaces is intrinsically interdisciplinary. This review will high- 
light computational techniques which are trivial to either the solid or liquid state community 
for the sake of bridging the gap between them. Important modeling studies on bare electrode 
surfaces have been made on crystal facet effects, surface reconstructions and terminations, 
non-stoichiometric compositions, and defects.—""— More effort in this area is needed to eluci- 
date surface structures which are starting points for interfacial studies. But the electrolyte 
is undeniably important and has been somewhat neglected. To take an example, DFT- 
computed energy dffererences for Li insertion into bulk cathode materials from Li metal 



have been reported as "open circuit voltages" (OCV),— which appears a misnomer ("intrin- 
sic potential" would be more appropriate). Measuring OCV invariably involves immersing 
an electrode into a liquid electrolyte. The absolute values of such voltages are modified 
by net surface charges, interfacial dipoles, and even purely quantum mechanical interfacial 
effects.—"— Even if interfacial voltage effects are small in magnitude for battery electrodes, 
they may become crucial for atomic level study of Li+ insertion and SEI formation. We 
will present other examples where explicit depiction of molecules qualitatively alters our 
understanding of battery degradation processes. Current theoretical studies on organic 
solvent-based Li-air batteries^^ have also focused on either the solid^~— or liquid^ state. 

Modeling reactions at LIB electrode/electrolyte interfaces is arguably a fledging, special- 
ized area. However, it is part of the broader fleld of computational and theoretical elec- 
trochemistry, intimately connected to aqueous interfaces,—"— notwithstanding the fact that 
water itself must be excluded from LIB. DFT-based computational studies of water- material 
interfaces pertinent to electrocatalysis,— ~— energy conversion (water splitting) ,—1^"— geo- 
chemistry and mineral dissolution,—"— and fuel cells^"— can inform and be informed by 
LIB simulations. Our work in fact owes part of its motivation^ to water dissociation studies 
on Ti02 surfaces.—"— In that example, the interactions of a sub-monolayer of water with 
different facets of anatase and rutile have been considered in joint DFT and experimental 
studies on clean surfaces under UHV conditions^"— before AIMD studies of water-Ti02 in- 
terfaces are conducted — including simulations with anions which may not be stable under 
unsolvated UHV conditions. We argue that a similar recipe may benefit fundamental studies 
on LIB. 

At the same time, it must be stressed that any frozen (i.e., UHV-like geometry opti- 
mization at T=0 K) description of electrolyte molecules is an uncontrolled approximation 
of liquids, although it can yield very useful insights. ^^ Theoretically, one usually claim that 
Af^aB if a low order perturbative expansion between A and B yields accurate predictions. 
Liquids and crystalline solids are separated by discontinuous first order phase transitions. 
Elementary calculus dictates that no real perturbative series connect them. In that sense, 
liquids and crystalline solids are not "similar;" they are profoundly different. As an example, 
salt solubility drops many orders of magnitude across the freezing point - yet salts are vital 
to electrochemistry! Solvent oxidation/reduction involves changes in solvent charge states, 
and the stabilization of these reaction products is related to ionic solubility. Modeling an 



electrode with a net charge is facihtated by exphcit treatment of hquid electrolytes that 
compensate the charge. AIMD with explicit molecular description of liquids is rigorously 
suited for modeling salt effects at interfaces, although this remains costly at present. 

This article is organized as follows. Sec. HT] and UTTl describe electrolyte decomposition re- 
actions on pristine graphite and lithium metal surfaces at the initial stages of SEl formation, 
respectively. We observe two-electron induced ethylene carbonate decomposition via break- 
ing two different C-0 bonds, instead of the one-electron pathway which has been the focus 
in the literature. Sec. IIVI discusses electrode coated with an insulator, where the long-range 
electron tunneling rate is found to depend on both the molecular species being reduced and 
its charge state. Sec. |V] reviews solvent decomposition modeling on spinel manganese oxide 
cathode surfaces and suggests that solvent breakdown and Mn dissolution may be related. 
Sec. I VI I summarizes the article. A supporting information (S.l.) document is included to 
discuss computational challenges to be overcome and possible new directions of research, 
including the prediction and control of voltages, electron transfer, and DFT accuracy is- 
sues. In each section, we briefly describe the motivations and main predictions at a model 
LIB interface, discuss the significance in some detail, and bring up pertinent computational 
issues. 

II. LIQUID EC ON LICe EDGES: SEI FORMATION 

LIB charging occurs at almost —3 V vs. standard hydrogen electrode (almost V vs. the 
Li''"/L(s) couple) with anodes like graphite and silicon. Most electrolyte molecules are elec- 
trochemically reduced at such voltages. On graphitic carbon, if EC (Fig. [T)d) is a main 
component of the electrolyte, stable passivating SEI films are formed from electrolyte de- 
composition products, and they prevent further electron leakage to the electrolyte which 
causes its continuous breakdown.- This is somewhat analogous to passivation of some metal 
surfaces by oxide or hydroxide in contact with water. Given its critical importance and 
ubiquity, it is tempting to call EC the "new water" of lithium ion batteries. 

Experimentally, the battery is assembled without pre-intercalation of Li into graphite. 
During first charge, as the voltage is lowered (effectively putting a negative charge on the 
anode compensated by cations at the interface), SEI starts to form. At even lower voltages, 
Li''" begins to insert. This dynamical process is difficult to simulate partly because the excess 



charge/voltage relation is difficult to control. Long trajectories are needed to equilibrate and 
converge even the open circuit voltage drop between the interior of the electrode and the 
electrolyte region outside the double layer™ such a calculation may be slightly beyond cur- 
rent AIMD time scale,— especially for the purpose of equilibrating salt diffusion. Modeling 
an applied voltage in a completely condensed-phase setting is even more challenging. See 
the S.I. for a more detailed discussion. 

Instead, we have chosen to start with fully lithiated LiCg with a charge-neutral simulation 
cell (apart from a positive charge arising from one Li"*" ion in the electrolyte region). While 
this LiCe stoichiometry does not truly reflect the experimental conditions during initial 
charging, it allows us to set the Li chemical potential of LiCg in the solid state anode to 
the solid state value by approximately matching the energy of the last Li atom added to 
the model electrode, in the absence of electrolytes, with the chemical potential of Li in bulk 
LiCg. Unlike on inert electrodes, the voltage in the anode interior can be tuned by varying 
Li content without inducing a net charge, although the anode surface may still retain a new 
charge; the potential of zero charge has not been determined experimentally or theoretically 
for most intercalation charge-neutral materials. Thus, with this approach, there remains 
some ambiguity in the voltage due to liquid-solid interfacial effects (see the S.I.). In the 
next section, we show that the electrochemical reactions and products discussed in this 
section are general; they are also observed at Li metal interfaces, where (unlike graphite) 
one does not have to worry whether the SEI forms before Li intercalation. 

Figure Ell depicts an AIMD snapshot at the initial stages of EC decomposition on graphite 
pre-intercalated with Li.iS 32 EC molecules are sandwiched between the anode surfaces. 
They are pre-equilibrated using non-reactive classical force fields. 7 ps into the AIMD tra- 
jectory, two electrons from the electrode have been transferred to each of three EC molecules 
coordinated to Li"*" ions in bulk solution or at the graphite edge. These "EC^~" species de- 
compose by breaking Ce-0 and Cc-0 bonds (Fig.[2|i) into two sets of products, respectively: 

EC + 2e- -^ col- + C2H4; (1) 

EC + 2e- -^ OC2H2O2" + CO. (2) 

Consistent with these reactions, significant amount of C2H4 and CO gases have been detected 
in gas chromotography measurement during battery charging,—"— although the precise gas 
speciation varies with experimental groups.—"^ Carbon labeling techniques have shown that 
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FIG. 2: (a) EC liquid with a solvated Li"*" ion confined between pristine C=0 graphite edges after 
a 7 ps AIMD trajectory. Intact EC molecules are shown as wireframes and graphite sheets are 
depicted as stick figures, (b) EC liquid confined between Li metal surfaces after a 10 ps trajectory. 
Eleven instances of CO-route decomposition and one of COg" route are observed, (c) EC~:Li+ 
with a broken Ce-0 bond. Adding a second e~ to this species yields C2H4 and COg^, while 
recombination of two ring-opened EC~ radicals can form butylene dicarbonate. (d) EC^^:Li"'" 
with a broken Cc-0 bond; this is the weaker (lower barrier) bond for EC with two excess e~ 
and its cleavage leads to the CO-route. (c) and (d) are optimized using a dielectric continuum 
approximation.— 

EC is a significant source of CO.^^ COg^ is known to be a SEI component. OC2H40^~ is 
reactive and can form other products, including oligomers which can be further reduced.—*^ 
These AIMD predictions are significant for the following reasons. In the literature, the 
slower one-electron, Ce-0 bond breaking decomposition route has been much quoted. It 
has been invoked to explain most observed SEI products from ethylene dicarbonate (EDC, 
reportedly a major SEI component^) to oligomers.— 1^ Modeling efforts, including DFT 
cluster calculations^i and reactive force-field construction,— have also focused on this l-e~ 
Ce-0 cleavage route. Even if 2-e~ processes are discussed, as they should be at the initial 
stages of SEI growth when e~ transfer is fast, most published works refer to mechanisms 
that involving breaking Ce-0 bonds via Eq. [H CO gas release (Eq. [2]), reported in several 
measurements^'"— and consistent with breaking Cc-0 bonds (Fig. |2]i), have arguably been 
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somewhat ignored. Because our model system contains a large electron reservoir and many 
EC molecules, it supports multiple EC reduction reactions and decomposition products 
anchored at the electrodes. Our unbiased AIMD simulations, which do not dictate a prior 
pathways, predict both 2-e~ induced reactions (Eq. [1] and |2]) at this out-of-equilibrium, 
highly driven initial stage of SEI formation. Thus Cc-0 bond-breaking (Eq. [21 Fig. [2li) is 
at least as fast as Ce-0 (Fig. ^) cleavage. This has been confirmed in a quantum chemistry 
cluster-based publication^^ which also suggests that l-e~ processes only dominate when the 
SEI thickens and e~ transfer significantly slows down. A central insight of AIMD and AIMD- 
inspired studies is that the Cc-0 bond becomes very weak after EC has absorbed one or 
two electrons. Therefore Eq. [2] needs to be considered when interpreting experiments and 
in force field constructions. We believe the 2-e^, CO-gas route is a more logical mechanism 
to yield EDC, reportedly a main SEI component, than l-e~ mechanisms.— 

We reiterate that the voltage drop between the model electrode and electrolyte regions is 
not precisely controlled in these simulations, unlike in experiments. However, the intrinsic 
reduction voltage of intact EC~ has been predicted to be less negative than that of charge 
neutral EC within the accuracy allowed by the dielectric continuum approximation used 
therein."2^ This suggests that the two-electron route and can occur at any applied voltage 
where EC~ is formed. The precise applied voltage may still alter the electron transfer rate 
via Marcus theory and indirectly change the SEI product distribution.— 

Computational aspects: AIMD simulations involve solving, in real time, Newton's second 
law of motion F = ma using forces generated by DFT. By the ergodic hypothesis, a suf- 
ficiently long trajectory visits all pertinent liquid state configurations and allows sufficient 
sampling of equilibrium properties. A 10 ps AIMD trajectory is at least 10-100 times more 
costly than zero temperature optimization calculations for models of the same size. But 
when a finite temperature, explicitly liquid-state, all-molecule description of electrolyte is 
used, MD is the rigorous method. AIMD allows the use of theoretical models with the 
highest fidelity to experimental conditions to date. 

These AIMD simulations are based on the approximate DFT/PBE functional^ which 
tends to underestimate reaction barriers. Despite this, they suggest that much faster 2- 
e~ reactions will be observed than previous studies of l-e~ processes have suggested.— 
Our simulations of EC are conducted at T=450 K. One reason is that EC is a solid at 
room temperature. In batteries, the presence of co-solvents like DMC (Fig. [Tfc) reduces 



the melting point and viscosity. DMC cis-trans isomerization is currently beyond AIMD 
time scale, but advanced sampling technique may circumvent this problem in the future. In 
general, AIMD may not reproduce the correct temperature scale for liquid state structures 
and dynamics.—'^ In the celebrated case of liquid water, AIMD/PBE simulations at an 
elevated T=400 K is needed to yield T=300 K experimental water structure.— 

III. EC ON LI METAL SURFACES: AIMD AND UHV MODELING 

EC decomposition products on Li metal surfaces and on LiCg are qualitatively similar, 
although product compositions and salt effects differ in quantitative ways.— We have also 
conducted AIMD simulations of the liquid EC/Li metal (100) interface (Fig. [2b)-— "^ Here 
EC reduction reactions are even more violent. Within 10 ps, all 12 EC molecules touching 
the Li metal have reacted. 11 of the 12 release CO molecules react via Eq. [2]and one EC has 
decomposed into C03~ and C2H4 (Eq. [1]). Despite the use of a thermostat in the simulation, 
the heat generated from the reactions has melted the small Li model electrode. 

One signficance of this study is that the open circuit voltage (OCV) should be unam- 
biguously that of Li"''/Li(s). In contrast, in Fig. [2k, the voltage is estimated using solid 
state approximations, neglecting interfacial effects. -"^^ Computing the OCV remains a major 
challenge of computational electrochemistry (see the S.I.),—"— although perturbative ap- 
proaches (e.g., adding voltage-induced changes in Fermi levels for metallic electrodes) have 
been successfully applied at T=0 K.^^ Despite this ambiguity associated with the graphite 
anode, the same SEI products are predicted (Fig. [2^ and[2)D), showing the robustness of 
previous AIMD predictions.— "^^ 

Second, the smaller system size of Li electrodes permits the use of a more accurate but far 
more computationally costly hybrid DFT functional, namely HSE06,— in AIMD simulations. 
While only sub-picosecond trajectories are possible with AIMD/HSE06, EC is predicted to 
break the Cc-0 bond as before,— giving us further confidence in AIMD/PBE simulations. 
We have also found that increased Brillouin sampling in this smaller system does not yield 
different EC reduction reactions. 

Finally, the crystalline Li metal (100) surface enables straight-forward UHV condition 
DFT/PBE studies of single molecule EC decomposition. Using T=0 K nudged elastic band 
(NEB) calculations, the COg" route is found to be more exothermic than the CO route on 
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this surface.— Both mechanisms are shown to be almost barrierless, consistent with obser- 
vation of both reactions in picosecond AIMD/PBE trajectories (Fig. Wp)- The preference 
for the CO-route is attributed to kinematic factors. NEB calculations show that electro- 
chemical EC decomposition mechanisms are predicted to be viable for single EC molecules 
at low temperature on Li(lOO). This suggests that future UHV imaging of sub-monolayer 
EC at low temperature will be valuable. (Imaging molecules at liquid/solid interfaces is 
far more difficult.) Note that XPS measurements have been performed on 6-10 nm thick 
DMC®^ and 10-20 nm thick propylene carbonate (PC)^ films on Li metal surfaces under 
UHV conditions. Carbonates and alkoxides are observed in DMC films, suggesting both Ce- 
O and Cc-0 bond breaking, consistent with our findings for EC. PC differs from EC only 



by a methyl group. Alkyl carbonate products are reported in Ref. |86|. In our simulations, 
neither CO gas nor OC2H40^~ is a final product. CO is absorbed into the Li metal, while 
the end groups of OC2ll40^~ are reactive.—*^ Future interpretations of XPS measurements 
of EC decomposition may benefit from consideration of Cc-0 bond cleavage. 

Computational aspects: The NEB method, widely used in solid state physics/materials 
science under UHV conditions, is similar in spirit to transition state calculations performed 
using the Gaussian suite of programs. The transition state obtained should exhibit only 
one imaginary frequency along the reaction direction. Unlike solids at T=0 K, liquid state 
atomic configurations are not at local minima, and the instantaneous vibrational frequency 
spectrum computed at any MD snapshot contains many imaginary modes.— Therefore NEB 
cannot be used in the presence of an explicit liquid component (see Sec. |V]). Counter ions 
like PF^ are more difficult to stabilize under UHV conditions, and are better left to AIMD 
interfacial simulations which explicitly provide dielectric solvation.— 

IV. ATOMIC LAYER DEPOSITION AND ELECTRON TRANSFER 

We next consider a LiCe electrode strip coated with a 7 or 10 A thick, hydroxyl- 
terminated /3-LiA102 oxide layer (Fig. |3]).— These models mimic electrodes coated by the 
atomic layer deposition (ALD) technique.— ALD has been used to generate conformal, sub- 
nanometer thick oxide layers on electrodes with sub-Angstrom precision. This "artificial 
SEI" approach has shown technological promise in limiting electrolyte decomposition and 
improving anode and cathode cyclability. But ALD oxide layers are also ideal surrogates for 
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(a) 

FIG. 3: (a) A single EC~ adsorbed on LiA102-coated LiCg under ultra-high vacuum conditions. 
Electron transfer estimates are made based on this modelJS (b) EC liquid in contact with LiA102- 
coated LiCg strip. Al is depicted as yellow. 

insulator films, including SEI naturally formed from electrolyte decomposition, in studies of 
e~ tunneling. 

We have estimated the "non-adiabatic" electron tunneling rate from the LiCe electrode 
to the adsorbed EC molecule at T=0 K using— "^ 
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where Vq is the e~ tunneling matrix element, A is the reorganization energy, AG° is —Qe^ 
added to the applied voltage, and $ is the reduction potential. We set AG=0, which mimics 
the very initial stage of SEI formation as the applied voltage is lowed, and estimate that 
A~2 eV using constrained density functional theory (cDFT, see below).— "2^ A is found to 
play a key role in limiting the e~ tunneling rate to about 1/s. A rigorous comparison of 
apparent electron transfer rates associated with rapid EC decomposition on uncoated anode 
surfaces (Fig. [2^ & b) should in principle be made at the same applied voltage. However, we 
note that C=0 and C-OH terminations of LiCe model electrodes, which may yield slightly 
different open circuit voltages, yield similar electron-tunneling-induced EC decomposition 
rates.— 

Ref. 120 appears the first theoretical work to use a Marcus theory-like formulation (Eq. [3]) 



to study electron transfer through insulating films on battery anodes. This work underscores 
the importance of depicting explicit molecules. It is incorrect to think of e~ tunneling purely 
as a quantum mechanical barrier-crossing problem, as though the electrolyte were a feature- 
less e~ sink. Instead, EC and other electrolyte molecules/ions exhibit distinct $ and A that 
may lead to vastly different reduction rates. Perhaps more significantly, on a per molecule 
basis, EC is predicted to take on a second electron (EC~— >EC^~) at a rate much faster than 
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the first reduction event (EC— t-EC"),^^ although much more rigorous apphcations of Marcus 
theory should be used to improve e~ tunneling rates in the future. From this insight and 
using cluster-based quantum chemistry calculations, we have estimated, speculatively but 
apparently for the first time, the crossover between l-e~ and 2-e~ processes.™ The SEI has 
often been described as consisting of an inner, inorganic layer and an outer layer made of 
organic carbonates.-""- In our picture, the inner layer should consist of 2-e~ products, namely 
COg" and compounds arising from subsequent reactions of the reactive OC2H40^~ (Eq. [1] 
and[2]). The outer layer is likely butylene dicarbonate arising from barrierless recombination 
of l-e~ induced ring-opened EC~ radicals (Fig. Efc),— but oligomers^ and products arising 
from proton transfer between EC~ and intact EC molecules may become competitive when 
e^ transfer to the electrolyte is slow. 

A key difference between SEI growth and aqueous electrochemical reduction processes 
is that, in LIB, e~ can be transferred to the solvent molecules (EC and DMC) right at 
the electrode surface. In water, e~-accepting ionic species are often well-solvated ("outer- 
shell") complexes located at least Angstroms away from the electrode, and they experience 
stronger screening of electric fields by electric double layers. However, other electrochemical 
properties computed using AIMD in aqueous electrolyes, such as surface potential^'S^ and 
reduction potentials,—!^ are important for organic solvents. In fact, predicted EC reduction 
potentials^ can deviate by tenths of eV from the oft-quoted 0.7 to 0.8 V onset observed 
in experiments on graphite and even on Ti02 anodes.— More accurate predictions of $,— 
especially for EC adsorbed on electrode surfaces, will improve the predicted e~ transfer rate 
via the AG term in Eq. |3l Using explicit solvent models to calculate $ should also be 
considered. 

Computational aspects: Electron transfer may be "adiabatic" or "non-adiabatic" (Eq. [3]) 
depending on Vo and and A.— "^ In the former regime, electronic configuration responds 
instantaneously to nuclear degrees of freedom. DFT calculations assume such a Born- 
Oppenheimer separation of time scales, and are appropriate at the initial stages when the 
electrolyte is in contact with pristine metallic electrodes. In the non-adiabatic regime, elec- 
tron tunneling is slow, and the system is not necessarily in its ground electronic state. DFT 
is problematic here. For example, it cannot confine e~ behind an insulating film if electrons 
are only metastable on the electrode at a particular atomic configuration. 

Thus ALD-coated electrodes present a challenging prototype problem for DFT-based 
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studies. We have found that T=0 K DFT NEB calculations based on the widely used PBE 
functional allows the unphysical splitting of an electron between the electrode and the EC 
molecule without a large energy cost.— This appears to be an example of self-interaction 
error (SIE).^^ Consequently, the e~ tunneling barrier appears underestimated using this 
functional. The problem may be alleviated using more costly hybrid functionals which are 
less susceptible to SIE. 

We have instead applied a combination of PBE calculations and cDFT which has been 
revitalized by van Voorhis and others.—"^ We have implemented local shape functions fi{r) 
centered around all atoms ( "i" ) of the EC molecule, which allows us to approximately add 
or subtract one e~ into or from a EC molecule adsorbed on the insulating LiA102 film 
at T=0 K. This gives A predictions via relaxation of vertical electronic excitations,— >^i22 
and yields an estimate of the non-adiabatic Vo as a by-product.— i^^S Preliminary attempts 
to compute vertical excitations at explicit liquid-solid interfaces (Fig. [3]d), related to finite 
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and the S.I. 



temperature A calculations, have been made but they remain costly. See Ref. 

for discussions of the approximations used and avenues of improvement in this important 

area. 

Long-range electron transfer is clearly pertinent to aqueous electrolyte interfaces where 
the electrode is coated with an oxide or hydroxide layer. So far most computational 
work on electron transfer at water-material interfaces have focused on pristine metal 
electrodes.— »^"-i2^ Many energy applications feature non-precious metal electrodes that re- 
act with the solvent to form insulating films. 

V. EC DECOMPOSITION ON CATHODE SURFACES 

We have also studied EC decomposition on the (100) surfaces of spinel Lia;Mn204,— a 
promising choice of LIB cathode material.- The issues to be addressed are twofold. First, 
electrolytes are known to decompose on spinel oxide surfaces despite the fact that the working 
voltage range of Lia;Mn204 is not high. The "intrinsic" oxidation voltage of EC molecules, 
computed in the absence of electrodes even under conditions most favorable to oxidation 
(EC coordinated at PFf^ in a low dielectric constant mediumr^)^ lies outside the spinel 
operational window. This suggests that, for EC to be oxidized, the surface has to play 
a catalytic/reactive role. The adsorbed, oxidized organic fragments have not been fully 
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FIG. 4: (a)-(b) Potentials-of-mean-force for two segments of the EC oxidation reaction, (c) Intact 
EC (configuration C) on Lio.6Mn204 (100) surface, (d) Intermediate D. Note that the EC Cc atom 
sits atop a surface oxygen ion that is not bonded to a Mn immediately below, (e) Intermedate E, 
with a broken Cc-Oe bond. The surface Mn(III) ion coordinated to the Oe now becomes a Mn(IV). 
(f) Product F; a proton and two electrons are transferred to the surface. 

identified. Second, the spinel oxide also degrades and Mn ions dissolve from it.- 

To our knowledge, Ref. l2l| is the first theoretical study of organic solvent molecule reac- 
tions on LIB cathodes. Unlike on anodes, unconstrained AIMD trajectories of EC liquid on 
defect-free Lio.6Mn204 (100) surfaces do not lead to spontaneous electrochemical reactions. 
The required timescale is apparently too long.^"^ To guide AIMD studies, we first consider 
a pristine crystal surface under UHV conditions. NEB calculations show that chemisorption 
of a single EC with an internal Cq-Oe bond broken and with an EC oxygen atom strongly 
bound to a 5-coordinated surface Mn(IV) ion preceeds oxidation. This chemisorbed geom- 
etry leads to a very exothermic transfer of two electrons and a proton to the oxide surface. 
DFT/PBEO-predicted energetics are similar for the intermediates.— 
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Then AIMD potential-of-inean-force (PMF, see below) simulations are conducted on ex- 
plicit liquid EC/electrode interfaces at finite temperature to investigate this UHV-motivated 
pathway. These simulations yield results similar to UHV predictions (Fig. Hj). In the first 
two steps, no charge transfer occurs, and the dielectric solvation missing in UHV studies 
should not play a strong role. In the key third step, electrons and a H+ are transferred, and 
the liquid environment slightly lowers the reaction barrier compared to UHV calculations. 
The predicted barrier appears readily surmountable in battery operation time scale. 

One significance of this work is that it suggests solvent oxidation and electrode degrada- 
tion may be related. The oxidized EC molecule fragment has pulled an oxygen ion out of the 
surface (Fig. Hf). The proton transfer from EC to the oxide is perhaps even more intrigu- 
ing. Acid has been found to accelerate Mn(II) dissolution.^ Acid protons have often been 
assumed to come from trace water, inevitably present) reacting with the PF^ salt. However, 
Oh and coworkers have reported evidence that protons can come from organic solvent, espe- 
cially tetrahydrofuran with LiC104 salt.-^ Our proposed EC oxidation mechanism may take 
place cooperatively with the trace water route, because H+ deposited on the oxide surface 
may create H2O as Mn dissolves. This work paves the way for future studies of oxidation 
on spinel (111) surfaces, which are more prominent than (100). Incidentally, the surface 
reconstruction of the (111) surface in vacuum has not been published.— For liquid-state ex- 
perts, this last point emphasizes the difficulty of identifying appropriate starting surfaces to 
perform liquid-solid interfacial calculations. Finally, as-synthesized LiMn204 particles are 
covered with lithium carbonate and inactive surface oxide phases. For Mn ions to dissolve, 
the carbonate film must first either crack-i2I or dissolve in the electrolyte. 

Much attention has been paid to Mn charge states, especially on bare spinel oxide 
surf aces. -^i^ This is because Mn(II) dissolution has been proposed to occur via two Mn(III) 
disproportionating into Mn(II) and Mn(IV).^ However, for an Mn ion to dissolve directly 
from the spinel surface, it must first coordinate to one or more solvent molecules.— We have 
observed that, when an EC molecule binds to a 5-coordinated surface Mn(III) on the (100) 
surface under UHV conditions, the now 6-coordinated Mn(III) transforms into a Mn(IV) by 
donating an electron to a subsurface Mn(IV). This seems to render Mn charge states on the 
bare (100) surface irrelevant, and emphasizes that molecules must be explicitly depicted to 
understand Mn dissolution. 

Computational Aspects: As discussed in Sec. IIIIl NEB cannot be used to find reaction 
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barriers in the presence of explicit liquid components. Instead, AIMD PMF calculations are 
performed using the umbrella sampling method.-i^S, Reactive coordinates R that continuously 
link the reactant, transition states, and product are chosen using UHV results as guide. The 
free energy change along R is proportional to —k^TlnPlR), where P{R) is the probability 
that the value R occurs in a trajectory after correcting for effects of umbrella constraints 
(Fig. H^ & b).- - Calculating the reaction free energy (AG), which is path-independent, also 
benefits from having a continuous R path. An alternative would be to calculate the average 
enthalpies of the initial and final state and subtracting them. This alternate approach is 
plagued by thermal uncertainties which scale as {NT/t)^'^kBT, where A^ is the number of 
nuclear degrees of freedom, t is the total trajectory length, and r is the correlation time. 
If the simulation is faithful to experiments and features many liquid molecules (large A^), 
the noise can easily be a fraction of an electron volt. The PMF method circumvents this 
problem. 

VI. CONCLUSIONS 

We have reviewed our modeling work on lithium ion battery electrode/electrolyte interfa- 
cial electrochemical reactions leading to SEI formation on anode surfaces and the reactions 
between electrolytes and spinel lithium manganese oxide cathode surfaces. AIMD simu- 
lations are valuable for predicting unbiased, fast, kinetically-determined multiple reaction 
mechanisms at the initial stages of SEI formation. Although costly, insights gained using 
AIMD can motivate new research directions which require less computationally intensive 
methods. Electron tunneling between electrode and electrolyte through an insulating layer 
will be an important area of research (see also the S.I.). The explicit depiction of molecules 
and electrode surfaces is shown to be important for both electron transfer and cathode 
degradation. The utility of ultra-high vacuum experiments on clean electrode surfaces cov- 
ered with sub-monolayers of electrolyte molecules is proposed. We have also highlighted 
the synergy between lithium ion battery modeling and theoretical studies of water-material 
interfaces.—"^— Exchange of knowledge between these disciplines will accelerate progress in 
modeling battery processes where more fundamental modeling studies are necessary. 
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